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Abstract — Based on data for 102 OB3 stars with known proper motions and ra- 
dial velocities, we have tested the distances derived by Megier et al. from interstel- 
lar Ca II spectral lines. The internal reconciliation of the distance scales using the 
first derivative of the angular velocity of Galactic rotation Q' and the external recon- 
ciliation with Humphreys's distance scale for OB associations refined by Mel'nik and 
Dambis show that the initial distances should be reduced by ~20%. Given this correc- 
tion, the heliocentric distances of these stars lie within the range 0.6-2.6 kpc. A kine- 
matic analysis of these stars at a fixed Galactocentric distance of the Sun, R = 8 kpc, 
has allowed the following parameters to be determined: (1) the solar peculiar veloc- 
ity components (U Q , V Q , W Q ) = (8.9,10.3,6.8) ± (0.6,1.0,0.4) km s -1 ; (2) the Galactic 
rotation parameters — —31.5 ± 0.9 km s" 1 kpc -1 , Q' = +4.49 ± 0.12 km s -1 kpc -2 , 
S7q = —1.05 ±0.38 km s -1 kpc -3 (the corresponding Oort constants are A = 17.9 ±0.5 km 
s -1 kpc -1 , B = —13.6 ± 1.0 km s -1 kpc -1 and the circular rotation velocity of the so- 
lar neighborhood is |Vo| = 252 ± 14 km s -1 ); (3) the spiral density wave parameters, 
namely: the perturbation amplitudes for the radial and azimuthal velocity components, 
respectively, fn = —12.5 ± 1.1 km s -1 and fg = 2.0 ± 1.6 km s _1 ; the pitch angle for 
the two-armed spiral pattern i = —5.3° ± 0.3°, with the wavelength of the spiral density 
wave at the solar distance being A = 2.3 ± 0.2 kpc; the Sun's phase in the spiral wave 
X0 = -91°±4°. 



INTRODUCTION 

Data on various objects are used to determine the Galaxy's kinematic parameters. These 
include the radial velocities of neutral and ionized hydrogen with the distances derived 
by the tangential point method (Burton 1971; Clemens 1985; Fich et al. 1989), Cepheids 
with the distance scale based on the period-luminosity relation, open star clusters and 
OB associations with photometric distances (Mishurov and Zenina 1999; Rastorguev et 
al. 1999; Dambis et al. 2001; Zabolotskikh et al. 2002; Bobylev et al. 2008; Mel'nik 
and Dambis 2009), and maser sources with trigonometric parallaxes (Reid et al. 2009; 
McMillan and Binney 2010; Bobylev and Bajkova 2010). 
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Young massive luminous stars with large heliocentric distances are important for the 
solution of our problem. However, among, for example, the O-type stars from the Hippar- 
cos (1997) catalog, only 12 have parallaxes differing significantly from zero and the most 
distant of them, 10 Lac, is at a distance of only ^540 pc from the Sun (Maiz-Apellaniz et 
al. 2008). The distances to a large number of OB stars determined spectroscopically from 
the broadening of interstellar Call, Nal, or KI absorption lines are of indubitable interest. 
The method for determining such distances has long been known (for a review, see Megier 
et al. 2005). However, only recently did Megier et al. (2005) and Megier et al. (2009) 
tie the "calcium" scale, respectively, to the Hipparcos (1997) trigonometric parallaxes. 
Megier et al. (2009) estimated the accuracy of an individual distance to OB stars to be 
~ 15%. For many of them, such highly accurate distance estimates have been obtained 
for the first time, given that their Hipparcos trigonometric parallaxes are not significant. 
The method is based on the assumption about a uniform distribution of ionized atoms in 
the Galactic plane. As these authors point out, the derived distances are applicable only 
at low elevations of stars above the Galactic plane (\z\ < 0.8 kpc). They also point to 
possible local inhomogeneities in Galactic longitude, especially in the region of the cluster 
Tr 16. 

The best-studied inhomogeneities are those in the immediate solar neighborhood as- 
sociated with the Local Bubble. Here, the variations in the number density of Call ions 
reach one order of magnitude (at a mean number density ricaii ~ 10" 9 cm -3 ) and the 
spatial sizes of the inhomogeneities are ~60 pc (Welsh et al. 2010). 

Our goal is to test the distance scale derived from interstellar Call spectral lines and, 
subsequently, to use it to investigate the kinematics of OB stars, in particular, to construct 
the Galactic rotation curve and to determine the spiral density wave parameters. To refine 
the "calcium" scale, we apply the method for internal reconciliation of the distance scales 
using the first derivative of the angular velocity of Galactic rotation Q' (Zabolotskikh et 
al. 2002) and the method for external reconciliation with Humphreys's (1978) distance 
scale for OB associations refined by Mel'nik and Dambis (2009). 

DATA 

We used data on 290 young OB3 stars whose distances were determined by Megier et 
al. (2009) from the equivalent widths of Call K and Call H lines by tying them to the 
trigonometric parallaxes of a revised version of the Hipparcos catalog (van Leeuwen 2007). 
We supplemented the sample with data from the CRVAD-2 compilation (Kharchenko et al. 
2007), with contains the radial velocities, proper motions, and photometric characteristics 
for ^55000 stars. 

For spectroscopic binaries, we checked them against the SB9 database (Pourbaix et al. 
2004) in order to refine their systemic radial velocities Vy. For a number of stars, we made 
significant changes to the CRVAD-2 radial velocities. There are also the most recent V 1 
determinations, for instance, for such stars of our sample as HIP 31978 (Cvetkovic et al. 
2010), HIP 32067 (Mahy et al. 2010), or HIP 78401 (Tango et al. 2009). 

As a result, we obtained a working sample of 258 Hipparcos OB3 stars with the 
distances from Megier et al. (2009), proper motions (van Leeuwen 2007), and radial 
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velocities. These stars have various luminosity classes. Note that ~ 20% of the sample 
are either known runaway stars (Tetzlaff et al. 2011) or candidates for runaway stars due 
to their large (> 40 km s -1 ) residual space velocities. 

THE METHOD 

The method used here to determine the kinematic parameters consists in minimizing a 
quadratic functional F: 

N N N 

min F = J2wi(VJ - Vjf + £ w>(V> - V>f + - V{f (1) 

3=1 3=1 3=1 

where N is the number of stars used; j is the current star number; V r is the radial velocity; 
V t = 4.74r/x/C0sfc and V b = 4.74r(x b are the proper motion velocity components in the Z 
and b directions, respectively, with the coefficient 4.74 being the quotient of the number of 
kilometers in an astronomical unit and the number of seconds in a tropical year; Vj , Vf , V b J 
are the measured components of the velocity field (data); w-j., wj,wl are the weight factors, 
provided that the following constraints derived from Bottlinger's formulas (Ogorodnikov 
1965) with an expansion of the angular velocity of Galactic rotation Q into a series to 
terms of the second order of smallness with respect of r/R and with allowance made for 
the influence of the spiral density wave hold: 

V r = —u & cos b cos(/ — lo) — Vq cos b sin(Z — l ) — w & sin b — (2) 
—R (R — Rq) sin(/ — Z ) cos bQ' Q — 
-0.5Ro(R - R ) 2 sin(/ - Z ) cos6fi' / + 
+ve sin(Z — l + 9) cos b — vr cos(/ — l + 9) cos b, 
Vi = u Q sin(7 — l ) — v Q cos(/ — / ) — 
— (R — R )(R cos(/ — l ) — r cos6)fi' — 
-(R- R ) 2 (Ro cos(/ -l )-r cos 6)0.5fig + rfi cos b+ 

+v e cos(/ - lo + 9) + v R sin(/ -l + 9), 
V b = u Q cos(/ — Iq) sin b + v Q sin(/ — l ) sin 6 — w Q cos b+ 
+R (R - i? ) sin(Z - Z ) sin 
+0.5R {R - R ) 2 sin(Z - Z ) sin 6fig- 
— tig sin(Z — Z + 0) sin 6 + v R cos(Z — Z + 9) sin 6, 

where r is the stars heliocentric distance; the stars proper motion components Hi cos b 
and Hb are in mas yr _1 and the radial velocity V r is, in km s _1 ; Uq,v q ,w & are the stellar 
group velocity components relative to the Sun taken with the opposite sign (the velocity 
u is directed toward the Galactic center; v is in the direction of Galactic rotation; and w 
is directed to the north Galactic pole); Ro is the galactocentric distance of the Sun; R is 
the galactocentric distance of the star; l is the direction to the kinematic center (to the 
Galactic center) — we included it as an unknown to reveal possible peculiarities of the 
"calcium" stellar distance scale. f2 is the angular velocity of rotation at the distance R ; 
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the parameters Q' and Qq are, respectively, the first and second derivatives of the angular 
velocity; the distance R is calculated from the expression 

R = (rcos bf -2 J R rcos6cos(/ - l ) + R 2 . (3) 

To take into account the influence of the spiral density wave, we used the simplest kine- 
matic model based on the linear theory of density waves by Lin et al. (1969), in which 
the potential perturbation is in the form of a traveling wave. Then, 

vr = f R cos x, (4) 
v e = f e sin x, 
X = m[cot(i) HR/Rq) -9] + Xe , 

where /r and fg are the perturbation amplitudes for the radial (directed toward the 
Galactic center in the arm) and azimuthal (directed along the Galactic rotation) velocities; 
% is the spiral pitch angle {% < for winding spirals); m is the number of arms, with m = 2 
taken here; 9 is the position angle of the star (measured in the direction of Galactic 
rotation); Xo is the phase angle of the Sun, measured here from the center of the Carina- 
Sagittarius spiral arm (R « 7 kpc), as was done by Rohlfs (1977). The parameter A is 
the distance (in the Galactocentric radial direction) between adjacent segments of spiral 
arms in the solar neighborhood (the wavelength of the spiral wave) — it is calculated from 
the relation 

tanW = 2^iV (5) 

The described method of allowance for the influence of the spiral density wave was applied 
by Mishurov and Zenina (1999) and Fernandez et al. (2001), where its detailed description 
can be found, and by Zabolotskikh et al. (2002). 

The weight factors in functional (1) are assigned according to the following expressions 
(for simplification, we omit the index j): 



w r = S /JS$ + al, 



Wl = (3 2 S /y/SZ + <7 



2 

VP 



where So denotes the dispersion averaged over all observations, which has the meaning of 
a "cosmic" dispersion taken to be 8 km s -1 ; (3 = a Vr /a Vl and 7 = a Vr /a Vb are the scale 
factors that we determined using data on open star clusters (Bobylev et al. 2007), /3 = 1 
and 7 = 2. The errors of the velocities V\ and V& are calculated from the formula 



^(VLH) = 4 - 74r 



\ 




The optimization problem (l)-(4) is solved for eleven unknown parameters 
u & ,v Q ,w & ,flQ,flQ,flQ,lo,fji,fg,i, and x& by the coordinate-wise descent method (the 
sought-for parameters are taken as the coordinates). 
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Table 1: A summary of fig and X/Rq determinations 



Reference 


Method 


Data 


n o (Ro = 8 kpc), 


X/R 








km s _1 kpc~ 2 




Clemens (1985) 


AV e 


CO+HI 




0.22 


Mel'nik et al. (2001) 


V(fl) + V r 


OB associations 




0.28 ±0.03 


Zabolotskikh et al. (2002) 


V(n) + V r 


Cepheids+OSCs 


-4.3 ±0.2 (*) 


0.33 ±0.04 


Zabolotskikh et al. (2002) 


V(n) + v r 


Supergiants 


-4.4 ±0.2 (*) 


0.36 ±0.05 


Bobylev et al. (2008) 


AV 9 


HI+HII+OSCs 




0.33 ±0.04 


Bobylev et al. (2008) 


Vr 


OSCs 




0.23 ±0.07 


Mel'nik and Dambis (2009) 


V(H) + V r 


OB associations 


-4.4 ±0.2 




Bobylev and Bajkova (2010) 


Vr 


Masers 


-4.5 ±0.2 


0.25 ±0.03 



* : The mean of the values obtained at Ro = 7.5 kpc and Rq = 8.5 kpc; we calculated X/Ro us 
from the pitch angle i (for the two-armed model) based on Eq. (5). 



We estimated the errors of the sought-for parameters through Monte Carlo simulations. 
The errors were estimated by performing 100 cycles of computations. For this number 
of cycles, the mean values of the solutions virtually coincide with the solutions obtained 
from the input data without any addition of measurement errors. Measurements errors 
were added to such input data as the radial velocities, proper motions, and distances. 

Here, we take a fixed value of Ro- Reid (1993) found a weighted mean from the mea- 
surements published over a 20-year period, R = 8.0 ± 0.5 kpc. Taking into account the 
main types of errors and correlations associated with the classes of measurements, Niki- 
forov (2003) derived the "best value", Ro = 7.9±0.2 kpc. A similar result was obtained by 
Avedisova (2005), Rq = 7.8 ±0.3 kpc. Note several most recent Ro determinations. These 
include the direct distance measurements based on the orbits of stars moving around a 
massive black hole at the Galactic center, which give R = 8.4 ±0.4 kpc (Ghez et al. 2008) 
or R = 8.33 ± 0.35 kpc (Gillessen et al. 2009). A summary of the latest determinations 
can be found in Foster and Cooper (2010), where the weighted mean is Ro = 8.0 ±0.4 kpc. 
Given all uncertainties, we consider R = 8.0 ± 0.4 kpc to be the most probable value. 

RESULTS 

When solving the system of equations (2), 38 of 258 stars are rejected according to the 3a 
criterion already at the first step. Seventeen of these stars are identified with the catalog 
of candidate runaway stars (Tetzlaff et al. 2011); most of the remaining ones are either 
new candidates for runaway stars or have measurements of the radial velocities or proper 
motions that are too unreliable. 

Since the Gould Belt has a significant influence in the region r < 0.8 kpc, we do not 
consider any stars from this region when studying the parameters of the Galactic rotation 
curve. 

There is only one star in the region r > 3.2 kpc, HIP 85020 (r = 4.5 kpc). It has a low 
residual velocity, but we do not consider this star, because it is too far from the common 
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Table 2: Kinematic parameters found using the refined distance scale {p sca ie = 0.8) 



JT dbL dlllc LCI a 


Nn 1 
IN U 1 


Nn 9 
IN (J Z 


Nn ^ 
IN U O 


Nn A 
IN U *± 


Nn ^ 
IN U O 


(Zq , Kill b 




Q Q _|_ n 7 

o.y it u. t 


y.o it u. / 


Qnj.ni; 
y.u it u.o 


Q O _|_ H fi 
y.Z it U.O 


c/0, Kill b 


1U.O It l.U 


iu.o it u.y 


in 1 -1-1 n 

1U.1 It l.U 


10 4 4-1 1 

1U.4: It J L 


10 9 4-1 
1U. Z It l.U 


1^0, km s ^ 


D.o it U.4 






D.O It U.Z 


a Q _|_ n A 
O.O It U.4 


On km s ^ Icnr — ^ 


-31 5 ± 9 

J- • KJ — 1 — KJ • \J 


—31 5 ± 8 

KJ J- • t_» — 1 — U • w 


—31.8 ± 1.0 


-29 8 ± 1 3 


—31 2 ± 9 

KJ J- • — 1 V J • \J 


S7 , km s _1 kpc~ 2 


4.49 ±0.12 


4.49 ±0.11 


4.53 ±0.13 


4.14 ±0.16 


4.46 ±0.12 


Qq, km s _1 kpc~ 3 


-1.05 ±0.38 


-1.05 ±0.36 


-1.15 ±0.42 


-0.93 ±0.34 


-1.12 ±0.40 


f R , km s" 1 


-12.5 ± 1.1 


-12.5 ±0.9 


-13.4 ±0.9 


-9.4 ±4.4 


-12.7 ± 1.1 


/ e , km s- 1 


2.0 ± 1.6 


2.1 ± 1.5 


0.7 ± 1.7 


2.5 ± 1.7 


1.8 ±1.6 


l , deg. 


-1.5 ± 1.3 


-1.5 ± 1.2 


-1.1 ± 1.3 


-2.3 ± 1.2 




i, deg. 


-5.3 ±0.3 


-5.3 ±0.3 


-5.4 ±0.3 


-5.3 ±44 


-5.4 ±0.3 


Xo, deg. 


-91 ±4 


-91 ±4 


-89 ±3 


-87 ± 47 


-91 ±4 


A, kpc 


2.3 ±0.2 


2.3 ±0.2 


2.4 ±0.2 


2.3 


2.4 ±0.2 


do, km s _1 


9.5 


9.5 


10.6 


8.3 


9.6 




102 


102 


102 


219 


102 



grouping of stars. As a result, our working sample consists of 102 stars. 
The Initial "Calcium" Distance Scale 

We obtained the following solution from 102 stars in the range of distances 0.8-3.2 kpc: 
(u , v Q , w q , ) = (9.3, 10.2, 8.4) ± (0.7, 1.3, 0.4) km s" 1 and 

fi = -28.6 ±0.8 km s" 1 , (6) 
n' Q = ±3.91 ±0.10 km s" 2 , 
nl = -0.81 ± 0.35 km s" 3 , 
f R = -11.9 ± 1.3 km s _1 , / e = 2±2kms- 1 , 
/ = -l°±2 o , i = -6.8° ±0.5°, x = -94°±5°, 

the error per unit weight is o"o = 11.2 km s^ 1 . The almost zero value of Iq shows that 
there are no significant global longitudinal peculiarities. 

Based on the pitch angle % found in solution (6) and Eq. (5), we find X/Rq = 0.37±0.03 
and then A = 3.0 ± 0.2 kpc. This value of A is too large compared to the determinations 
of this parameter by other authors using various independent data. This suggests that 
the "calcium" distance scale is stretched. To determine how much it is stretched, we use 
several kinematic methods for comparing the distance scales. 

Refining the "Calcium" Distance Scale 

We know a method (Dambis et al. 2001; Zabolotskikh et al. 2002) where the distance scale 
coefficient (p SC aie) is included as an additional unknown in the initial kinematic equations 
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012345 0123456 
rfCall),kpc 0.8*r (Ca II), kpc 

Figure 1: (a) Stellar distances 0.8 x r Humphreys versus distances r Ca ii', (b) stellar distances 
0.8 x r Humphreys versus distances 0.8 x r Ca ii', the dotted lines mark the outer boundaries 
of the sample of stars used to determine the kinematic parameters. 



(e.g., (2)). In this case, it can be determined by simultaneously solving the system of 
equations or can be found by minimizing the residuals according to the x 2 test. 

Unfortunately, this method did not give a reliable result in our case. The minimum 
of x 2 is reached at p sca ie ~ 0.3 — 0.4. The new distances are calculated as r NEW = 
Pscaie xr - However, for such a radical reduction of the initial "calcium" scale, the remaining 
model parameters lose any physical meaning. We assume that this is determined by the 
"calcium" scale calibration method — the nonlinear (hyperbolic) relation between the 
equivalent widths of Call spectral lines and trigonometric parallaxes (Megier et al. 2009). 
Thus, different values of the linear coefficient p sca ie should be applied for different ranges 
of distances. Other distance scale reconciliation methods are also known. 

Reconciliation of the derivative Q' . 

The method for internal reconciliation of the distance scales using the first derivative 
of the angular velocity of Galactic rotation Q (Zabolotskikh et al. 2002) consists in 
comparing the values obtained by separately solving the system of equations (2). For this 
purpose, we solved a simplified system of equations with the unknowns u Q ,v Q ,w Q and 
Qo,Q' ,Qq; when using the radial velocities, we fixed the parameters w Q and fl found 
from the proper motions. We found Q' Q (Vr) = 3.87 ± 0.35 km s _1 kpc" 2 only from the 
radial velocities and Q (V^) = 3.29 ± 0.26 km s _1 kpc~ 2 only from the proper motions. 
Then, p scale = 3.29/3.87 = 0.85. 

At present, there are quite satisfactory estimates of Q and A for external scale calibra- 
tion. A number of such results are presented in Table 1; they were obtained in different 
distance scales: the tangential point method (Clemens 1985), the Cepheid scale, the pho- 
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Figure 2: Positions of 219 OB3 stars in projection onto the Galactic plane. 



tometric scale (open star clusters — OSCs) reconciled with the Hipparcos trigonometric 
parallaxes, and the trigonometric parallaxes (masers). By comparing Q = 3.9 km s _1 
kpc~ 2 found in solution (6) with the mean from Table 2, we obtain p sca ie = 3.9/4.4 = 0.89. 

Reconciliation of the wavelength A. 

By comparing the mean from Table 1, X/Rq = 0.29, with the result of solution (6), X/Rq = 
0.37, we find p sca i e — 0.29/0.37 = 0.78. Having separately considered the Galactocentric 
radial velocities of our OB3 stars, we found X/Rq = 0.36 ± 0.03 (A = 2.9 ± 0.2 kpc). For 
this purpose, we applied the method of Fourier analysis described in detail previously 
(Bobylev and Bajkova 2010). In this case, p sca ie = 0.29/0.36 = 0.80. Thus, this method 
gives p sca i e close to 0.8. 

Reconciliation with Humphreys's refined distance scale. 

The membership in OB associations with the distances to them determined by Humphreys 
(1978) is specified in the catalog by Megier et al. (2009) for a considerable number of stars. 
A number of researchers of Humphreys's scale (Dambis et al. 2001; Mel'nik and Dambis, 
2009) conclude that it should be reduced by ~ 20%, i.e., the refined distances should be 
calculated as r NEW = 0.8 x r Humphreys . Figure 1 shows the "0.8 x r Humphreys - r Ca ii" and 
"0.8 x rHumphreys — 0.8 x fcaii" relations. In fact, the distance scale by Mel'nik and Dambis 
(2009) is along the vertical axis. It is clearly seen from the figure that the "calcium" scale 
should be reduced with a coefficient close to 0.8. 

As can be seen from Fig. lb, the stars with distances greater than 3 kpc already deviate 
significantly from the general relation even in the corrected "calcium" scale. We did not 
specially determined the outer boundary indicated in the figure, because it was revealed 
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Figure 3: Galactic rotation curve (solid line). The vertical line marks the position of the 
solar circle. The dashed lines indicate the la confidence intervals. 



automatically — the stars outside this boundary were either rejected according to the 3cr 
criterion when solving the system of equations (2) or had no velocity measurements. 



Galactic Rotation Curve 

Using several methods of analysis, we showed that the scale coefficient p sca ie nes within the 
range 0.7-0.9, with a mean close to 0.8. We use this mean to form the refined distances 
tnew — 0.8 x rcaii with which all of the subsequent computations are performed. 

The positions of 219 stars in the Galactic XY plane computed in the refined distance 
scale are shown in Fig. 2. 

Based on the sample of 102 relatively distant stars located in the range of distances 
tnew from 0.6 to 2.6 kpc, we found several solutions of the system of equations (2) 
given in Table 2. Solution 1 was obtained in the same way as solution (6) using all tree 
components V r , Vj, V&. To obtain solution 2, we used the same components V r , VJ, Vb but 
fixed the velocity w Q = 7 km s -1 . We found solution 3 only from two components, V r 
and Vi, with the fixed velocity w & = 7 km s -1 . Solution 4 was obtained in the same way 
as solution 1 using all three components V r ,V h V b from all 219 stars. 

Solution 4 shows that the influence of nearby stars degrades considerably the accuracy 
of determining the spiral-structure parameters. Solutions 1 and 2 are virtually identical. 
Comparison of solution 3 with solutions 1 and 2 shows that using the three components 
V r ,Vi,Vb is more preferable, i.e., the stars are not too far from the Sun for valuable 
information to be lost in the components V&. 

Solution 5 was obtained in the same way as solution 1 using all three components 
V r , Vi, Vb, but the direction to the Galactic center was fixed at Iq = 0°. All of the deter- 
mined parameters in solution 5 agree well with the results of solutions 1 and 2. Therefore, 
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Figure 4: Galactocentric radial velocities for 102 relatively distant OB3 stars (a), the 
sine wave with a period of 2.3 kpc and an amplitude of 12.5 km s" 1 (solid line), and the 
residual azimuthal velocities (b). The vertical line marks the position of the solar circle. 



below we will use the result of the more general case (solution 1). 

Figure 3 displays the Galactic rotation curve constructed according to solution 1 from 
Table 2. The dashed lines indicate the boundaries of the confidence intervals correspond- 
ing to a lcr error level calculated by taking into account the contributions from the error 
in the angular velocity and the uncertainty in R$. 

Figure 4 presents the Galactocentric radial velocities Vr (the direction away from the 
Galactic center is considered the positive one) for 102 relatively distant OB3 stars and 
the residual azimuthal velocities AV rot . The azimuthal velocities are residual, because the 
Galactic rotation curve found was excluded from them; both velocities Vr and AV ro t were 
freed from the group velocity. The sine wave associated with the influence of the spiral 
density wave was fitted into the radial velocities (based on solution 1 from Table 2). 

The Oort constants calculated using solution 1 from Table 2 are A = 17.9 ± 0.5 km 
s _1 kpc -1 and B = —13.6 ± 1.0 km s _1 kpc -1 . The circular rotation velocity of the solar 
neighborhood is |Vo| = 252 ± 14 km s _1 (calculated by taking into account the error in 
R of 0.4 kpc). 

Table 3 lists the stars with residual space velocities \V pec \ > 40 km s -1 . We calculated 
the velocities V pec by excluding the rotation curve found and the solar peculiar velocity as 
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well as by excluding the wave found in the radial velocities. These stars are either already 
known runaway stars or suitable candidates for runaway stars. 

DISCUSSION 

(1) The phase x& we found is almost equal to —n/2. Its value does not depend on the 
distance scale. The specific values of the solar peculiar velocity relative to the local 
standard of rest (LSR) depend on the phase Xq ( see Eqs. (19), (20) in Bobylev and 
Bajkova 2010). Since young stars experience perturbations from the spiral density wave, 
the components of their mean motion u Q and v Q can differ significantly from the velocities 
found from older stars (Dehnen and Binney 1998). Denote the solar peculiar velocity 
components unperturbed by the spiral wave by (Uq, Vq, W q )lsr- As is clearly seen from 
Fig. 4a, the sine wave passes through zero at it! = 8 kpc. Therefore, despite its significant 
amplitude (/#), the radial component of the wave causes no shift of the mean velocities 
(its main contribution to the velocity Uqlsr)- The amplitude of the azimuthal component 
of the spiral wave (fg) is almost zero (its main contribution to the velocity V &LSR ). Then, 
(Uq,V q , W q ) lsr = (u Q ,v Q ,w Q ) = (8.9,10.3,6.8) ± (0.6,1.0,0.4) km s" 1 according to 
solution 1 from Table 2. 

The present-day situation with the determination of this velocity by various methods 
but without invoking any data on young stars is as follows. Schonrich et al. (2010) took 
into account the stellar metallicity gradient in the Galactic disk and found the following 
components: (U e , V Q , W e ) LSR = (11.1, 12.2, 7.3) ± (0.7, 0.5, 0.4) km s -1 . Having analyzed 
the eccentricities of stars in the solar neighborhood, Francis and Anderson (2009) found 
(Uq, V q , Wo) lsr = (7.5, 13.5.2, 6.8) ± (1.0, 0.3, 0.1) km s _1 . Based on an updated version 
of the Geneva-Copenhagen survey (Holmberg et al. 2007), Koval' et al. (2009) found 
(Uq, Vq, Wq)lsr = (5.1, 7.9, 7.7) ± (0.4, 0.5, 0.2) km s -1 from stars born at a circumsolar 
distance by taking into account the radial migration of stars and the metallicity gradient 
in the Galactic disk. Thus, the solar peculiar velocity components are (Uq, Vq, W q )lsr — 
(10, 11, 7) km s -1 . The components we found from OB3 stars are in good agreement with 
these values. 

This is confirmed by the analysis of a sample of blue supergiants that are closest in 
evolutionary status to the OB3 stars of our sample: (Uq, Vq, Wo) lsr — (6, 11, 7) km s _1 
(Zabolotskikh et al. 2002), where the phase was x& = — 97° at a significant amplitude 
Jr = —6.6 ± 2.5 km s -1 . A similar result was obtained by Fernandez et al. (2001) both 
using a sample of OB stars, ([/©, V Q , W Q ) LSR = (8.8, 12.4, 8.4) ± (0.7, 1.0, 0.5) km s _1 , and 
from Cepheids, (U e , V e , Wq) LS r = (6.5,10.4,5.7) ± (1.2,1.9,0.7) km s" 1 , with a phase 
Xq ~ —30°, although the spiral wave amplitudes turned out to be insignificant. 

(2) The ratio of the spiral wave amplitudes we found from our sample of 102 OB3 
stars, /r = — 13 ± 1 km s _1 and fg = ± 1 km s _1 , is in good agreement with the 
analysis of blue supergiants by Zabolotskikh et al. (2002), f R = —6.6 ± 2.5 km s^ 1 and 
fg = 0.4 ± 2.3 km s" 1 , and a sample of Cepheids with close values of these parameters. 

At the same time, it can be seen from Fig. 4b that a wave with a significant amplitude 
(?»10 km s _1 ) but with a wavelength A ~ 1 kpc is present in the azimuthal velocities. As 
we see from the figure, the local standard of rest turns out to be shifted in the region 
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Table 3: Candidates for runaway stars 
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Note. R denotes a runaway star according to Tetzlaff et al. (2011). 
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R — > Ro in the direction opposite to the Galactic rotation. When a sample of masers 
was analyzed (Bobylev and Bajkova 2010; Stepanishchev and Bobylev 2011), a special 
allowance for the perturbations from the spiral wave in the region R — > Rq made it 
possible to estimate the velocity components (U Q , V Q , W q )lsr- 

(3) The parameters of the Galactic rotation curve Qq and Q' Q calculated from 102 
OB3 stars using the refined "calcium" distance scale (solutions 1,2,3, and 5 from Table 2) 
are in good agreement with the results of the analysis of Cepheids and blue supergiants 
(Zabolotskikh et al. 2002), OB associations (Melnik and Dambis 2009), and masers 
(Bobylev and Bajkova 2010; Stepanishchev and Bobylev 2011). 

CONCLUSIONS 

We tested the distances derived from the equivalent widths of interstellar Call spectral 
lines by Megier et al. (2009). For this purpose, we used a sample of 102 relatively 
distantOB3 stars with known proper motions and radial velocities. 

The internal reconciliation of the distance scales using the first derivative of the angular 
velocity of Galactic rotation Q' Q and the external reconciliation with Humphreys's distance 
scale for OB associations refined by Mel'nik and Dambis (2009), we showed that the initial 
distances should be reduced by ^20%. 

In the refined distance scale, our OB3 stars are located at heliocentric distances in the 
range from 0.6 to 2.6 kpc. We used them to construct the Galactic rotation curve and to 
determine the parameters of the spiral density wave. 

For a fixed distance to the Galactic center, Rq = 8 kpc, we found the solar peculiar 
velocity components (u e ,v Q ,w Q ) = (8.9,10.3,6.8) ± (0.6,1.0,0.4) km s -1 ; the angular 
velocity of Galactic rotation f2 = —31.5 ± 0.9 km s _1 kpc -1 and its derivatives fi = 
+4.49 ± 0.12 km s _1 kpc -2 , f^' = -1.05 ± 0.38 km s _1 kpc -3 . The corresponding Oort 
constants are A = 17.9±0.5 km s _1 kpc -1 and B = — 13.6± 1.0 km s _1 kpc -1 ; the circular 
rotation velocity of the solar neighborhood is |Vo| = 252 ± 14 km s -1 ; the amplitudes of 
the spiral density wave are f'n = —12.5 ±1.1 km s _1 and fg = 2.0 ±1.6 km s -1 ; the pitch 
angle of the two-armed spiral pattern is % — —5.3° ± 0.3° and the phase of the Sun in the 
spiral wave is x© = — 91°±4°; and the direction to the Galactic center is Zo = — 1.5°±1.3°. 
The wavelength of the spiral density wave at the solar distance is A = 2.3 ± 0.2 kpc. It is 
particularly pronounced in the Galactocentric radial velocities Vr. The residual azimuthal 
velocities AVg have a more complex structure. 
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